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We use the combination of configuration interaction and many-body perturbation theory method 
(CI+MBPT) to perform ab initio calculations the low-energy spectra of Crll with high accuracy. 
It is found that second-order MBPT diagrams should be included in a consistent and complete way 
for the MBPT to improve the accuracy of calculations in this five-valence-electron system. This 
contrasts with previous ions with fewer valence electrons where it was found that single-valence- 
electron diagrams dominate the corrections. Isotope shifts and relativistic shifts (g-values) are 
calculated for use in astronomical determination of the fine-structure constant in quasar absorption 
spectra. 
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I. INTRODUCTION 

Quasar absorption systems provide a unique probe of 
the value of fundamental constants throughout much of 
the visible Universe. The many-multiplet (MM) method 
enables the most complete analysis of optical spectra in 
the search for space-time variation of the fine-structure 
constant, a = e^/?ic [J, It makes use of all transi- 
tions seen in all ions in a given quasar absorption system 
to gain statistical significance and control systematics. 
Early results using spectra taken from the Keck telesco 
suggested that a may have been smaller in the past 

however when combined with new systems observed 
with the Very Large Telescope (VLT) the data is more 
consistent with a spatial variation in the fine-structure 
constant 6]. The gradient in values of a reconciles all 
existing measurements of a- variation Q- In particular 
the early Keck results that indicated a constant offset or 
"monopole" model, are entirely consistent with the spa- 
tial gradient "dipole" model since Keck mainly sees in 
the northern hemisphere (the a-dipole axis is oriented 
~ 30° from the equatorial axis). By contrast the VLT 
data is taken mainly in the southern sky. 

A spatial variation of a would manifest itself in a vari- 
ety of terrestrial 8] and astrophysical 9] systems, which 
could be used to confirm the dipole. It is also possible to 
devise complementary tests using subsets of the quasar 
absorption system data which may involve different sys- 
tematics. One such test, currently underway, is a variant 
of the many-multiplet method that only uses transitions 
in CrII and Znll [T^]. The transitions have opposite 
a-sensitivities and so a comparison of them is very sensi- 
tive to a- variation: Zn II transitions are s — p and hence 
their frequency increases if a increases, while the CrII 
transitions are d — p so their frequency decreases with 
increasing values of a. Furthermore the transitions are 
very close in energy. This means that only a small part 
of the optical spectrum is analysed, resulting in differ- 
ent (perhaps smaller) systematics. Of particular concern 
are "intra-order shifts" : velocity shifts of unknown ori- 
gin within each echelle order in the spectrograph [111, [13 ■ 



This systematic may differently affect measurements of 
a- variation when only CrII and Znll lines are utilized, 
compared to studies where a larger number and wider 
variety of transitions are used. 

One problem with using Cr II and Zn II transitions ex- 
clusively is that they are weak. Of course, this is the rea- 
son why they don't play a major role in the full MM anal- 
ysis despite being included whenever available. However 
there exist certain quasar absorption systems in which 
Znll lines are particularly strong [Tsj, and from these 
"metal strong" systems can be drawn a relatively large 
sample with which to perform the CrII/ZnII analysis. 

One potential systematic that has plagued all MM 
analyses is isotope abundance (l3 - [l^ . Isotopic structure 
cannot be resolved in the absorption spectra, so generally 
terrestrial isotopic abundances are assumed for the ab- 
sorber. Any deviation from terrestrial abundances would 
shift the centroid of the line profile, and this might mimic 
a change in a. Even quantifying the systematic can be 
difficult because the isotopic structures themselves are 
unknown for many of the UV transitions used in the MM 
analysis. The systematic is lessened in the context of a 
dipole result, since the isotope abundances would need 
to vary according to direction in the sky, which in itself 
would violate cosmological isotropy. Nevertheless, in or- 
der to quantify possible systematics the isotope structure 
should be known for all transitions used in the analysis, 
hence considerable efforts by many groups to calculate 
and measure them (see, e.g. |17h21| ). 

In this paper we calculate the isotope shifts and rela- 
tivistic shifts of the CrII transitions seen in quasar ab- 
sorption spectra. The corresponding parameters for Zn II 
have been calculated previously 15, 22,]. Our final results 
are presented in Tables |V] and IVII 



II. METHOD 

The ah initio CI-t-MBPT method 23] is described in 
full elsewhere ^3] ■ Details of relevance for our Cr II cal- 
culation are presented below. 
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A. Energy calculation 



Any perturbative theory works best when the pertur- 
bations are as small as possible. In CrII the d-wave 
electrons play an important role in shaping the atomic 
core, and so they should be included in the initial ap- 
proximation. As in previous works, our single-particle 
wavefunctions are calculated using Dirac-Fock (relativis- 
tic Hartree-Fock) . We explore two Dirac-Fock configura- 
tions: , which includes a half- filled 3cf subshell; and 
V^~^, which includes 3d^. In both cases we simply scale 
the potential due to the filled 3d subshell by the number 
of electrons to provide a "configuration averaged" ini- 
tial wavefunction. The choice of starting approximation 
is essential to obtaining a good final spectrum for CrII, 
but also leads to potentially large subtraction diagrams in 
many-body perturbation theory, as will be demonstrated. 

Once we have a Dirac-Fock potential for the core, we 
diagonalize the Dirac-Fock Hamiltonian over a set of 40 
B-splines {25l| spanning 40 atomic units to obtain a large 
set of valence and virtual orbitals from which we select 
those with the lowest eigenvalues. A set of configura- 
tions of valence electrons |/) are generated, from which 
eigenfunctions of the complete Coulomb-potential Hamil- 
tonian are calculated. We find that almost complete con- 
vergence of the CI calculation can be obtained using the 
basis 20spdf: that is we use s-wave states labeled 1 - 
20, p-wave states labeled 2 - 20, etc. (For the lowest 
eigenvalue states the label is just the principal quantum 
number.) With 5 valence electrons it is not possible to 
include all configurations and we must select those that 
contribute most to the wavefunction. We include all con- 
figurations that can be formed by one-particle excitations 
from the leading configurations 3d^, 3(i'*4s, and 3d''4p, as 
well as two-particle excitations from these same configu- 
rations up to the 6sp9d6f orbitals. The effects of higher 
orbitals and three-particle excitations were found to be 
small and were not included. 

Having achieved high saturation of the CI calcula- 
tion, core-valence effects are included using second-order 
MBPT by modifying matrix elements of the Hamiltonian: 



M 



{I\H\M) {M\H\J) 
E — Em 



(1) 



where the states \M) include all Slater determinants that 
have core excitations. The MBPT sum may be fur- 
ther separated into one-, two-, and three- valence-electron 
parts, denoted E^^), S^^) and E^^) in Refs. [Ulii]. Gold- 
stone diagrams and analytical expressions for these are 
presented in [23|. The states \M) include excitations from 
all core states into virtual states up to 30spdfg. The ef- 
fects on the energy calculation of including E^^^ E*^^\ 
and E*^'^^ sequentially are shown in Tables HI and HIl (for the 
V^~^ and starting approximations, respectively). 



TABLE I. Energy spectrum of CrII with orbitals calculated 
in the V^~^ approximation, relative to the experimentally 
determined ground state, ^S'5/2. Successive additions of 
E*^^ E'^', and E*^' are shown. (Energies in cm~\) 



Level 




J 


CI 


+E(^) 


+E(2) 


+E(3) 


Expt. [271 


3d*4s 




5/2 


8505 


-7294 


6888 


9491 


12148 






7/2 


8682 


-7128 


7095 


9699 


12304 




6po 


1/2 


41720 


27542 


40933 


44533 


46823 






3/2 


41805 


27635 


41028 


44638 


46905 






5/2 


41945 


27788 


41184 


44798 


47040 






7/2 


42140 


28000 


41401 


45036 


47227 






9/2 


42387 


28271 


41676 


45320 


47465 






11/2 


42687 


28600 


42008 


45668 


47752 




Gpo 


3/2 


43214 


29890 


42534 


46074 


48399 






5/2 


43333 


29982 


42686 


46216 


48491 






7/2 


43515 


30132 


42915 


46420 


48632 


Sd'^ip 


4po 


1/2 


43793 


30108 


42725 


46584 


48749 






3/2 


44057 


30328 


43038 


46907 


49006 






5/2 


44787 


30674 


43880 


47717 


49706 


Sd^-ip 




1/2 


44594 


29189 


43659 


47460 


49493 






3/2 


44664 


29323 


43741 


47552 


49565 






5/2 


44451 


29510 


43492 


47330 


49352 






7/2 


44840 


29721 


43959 


47727 


49646 






9/2 


45054 


29940 


44209 


47994 


49838 



TABLE II. Energy spectrum of Cr II with orbitals calculated 
in the approximation, relative to the experimentally de- 
termined ground state, 3d^ ^Ss /2 • Successive additions of E^^' , 



E'^', and E 


are 


shown. (Energ 


ies in cm 


-^) 




Level 




J 


CI 


+E(i' 


+E(2' 


+E(3' 


Expt. [27] 


3d*4s 




5/2 


6688 


-8662 


3953 


6415 


12148 






7/2 


6862 


-8515 


4157 


6617 


12304 


3d''4p 


6po 


1/2 


41540 


26023 


38794 


41627 


46823 






3/2 


41626 


26110 


38888 


41726 


46905 






5/2 


41768 


26257 


39044 


41887 


47040 






7/2 


41966 


26461 


39260 


42114 


47227 






9/2 


42217 


26724 


39534 


42399 


47465 






11/2 


42521 


27047 


39865 


42747 


47752 


3d'*4p 


Gpo 


3/2 


43124 


27918 


40888 


43534 


48399 






5/2 


43242 


27987 


41041 


43674 


48491 






7/2 


43421 


28113 


41265 


43885 


48632 


3d''4p 


4po 


1/2 


43816 


28030 


41092 


43995 


48749 






3/2 


44033 


28307 


41320 


44246 


49006 






5/2 


44812 


29073 


42251 


45175 


49706 


3d''4p 


6p,o 


1/2 


44432 


28788 


41749 


44708 


49493 






3/2 


44552 


28878 


41915 


44861 


49565 






5/2 


44303 


28603 


41591 


44537 


49352 






7/2 


44543 


28849 


41851 


44815 


49646 






9/2 


44766 


29041 


42108 


45078 


49838 



B. Isotope shift and relativistic shift 

Isotope shifts in atomic transition frequencies come 
from two sources: the finite size of the nucleus (field 
shift), and the recoil of the nucleus (mass shift). This 
mass shift is usually divided into the normal mass shift 
(NMS), which is easily calculated from the transition fre- 
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quency, and the specific mass shift (SMS). The mass shift 
is more important for hght elements, while for heavy el- 
ements the field shift dominates. In the case of CrII, 
the field shift is small; this paper is concerned with the 
mass-shift contribution, which is more difficult to cal- 
culate. The difference in the transition frequency, w, be- 
tween an isotope with mass number A' and an isotope A, 
Su}^ = oj^ — oj^^ can be expressed as [15] 



5lO ' — (^NMS + ^SMs) 



A', A 



(2) 



where Slr^^ is the change in mean-square nuclear charge 
radius. The normal mass shift constant can be expressed 
(in atomic units fi — e — = 1) 



1 X - 2 ^ 



where to„ = 1823 is the ratio of the atomic mass unit 
to the electron mass, and the sum is over all electron 
momenta, pj. The specific-mass-shift constant 



"-SMS — 

1<J 



2^ • Pj 



and field-shift constant F are more difficult to calculate. 
We use the non-relativistic form of the mass-shift opera- 
tor; relativistic corrections for optical transitions in light 
atoms are on the order of few percent and can be ne- 
glected [ll|. 

To calculate /csms we use the all-order finite-field scal- 
ing method. Here a rescaled two-body SMS operator is 
added to the Coulomb potential everywhere that it ap- 
pears in an energy calculation: 



1 



|ri - i'2| 



Api • P2 



We recover the specific-mass-shift constant as 



(3) 



(4) 



A=0 



The operator Q has the same symmetry and structure as 
the Coulomb operator (see Appendix A of Ref. We 
have previously shown that good agreement with exper- 
imental isotope shift can be obtained in many-valence- 
electron atoms and ions by u sing this finite-field in a 
CI-hMBPT energy calculation [11 [H, [H [H . 

The relativistic shift of a transition may be calculated 
in a similar fashion. We simply recalculate the transition 
energies, w, from the very beginning using different values 
of a from the laboratory value ao. The sensitivity to 
variation of the fine-structure constant is then extracted 
using 



duj 
da^ 



(5) 



TABLE III. Isotope shift constants fenMs and fcsMs for transi- 
tions to the ground state ^3^/2 (GHz-amu). 



Level 




J 


' " IN iVl O 


IvIS 

yN 


(CI) 

yN-l 


rvi s 

yN 


(CI-I-E) 

yN-l 


3d*4s 




5/2 


—200 


4520 


4326 


4944 


5062 






7/2 


—202 


4532 


4337 


4950 


5069 






1 /2 

1/ z 


—770 


3964 


3398 


4303 


4127 






3/2 


—771 


3970 


3403 


4314 


4140 






5/2 


—774 


3980 


3412 


4330 


4161 






7/2 


—777 


3992 


3423 


4346 


4178 






9/2 


-781 


4009 


3438 


4359 


4191 






11/2 


—785 


4028 


3456 


4362 


4188 




6po 


3/2 


—796 


4072 


3555 


4189 


4161 






5/2 


—797 


4080 


3563 


4143 


4090 






7/2 


-800 


4093 


3575 


4203 


4180 




Apo 


1/2 


-802 


3981 


3495 


4227 


4068 






3/2 


-806 


3983 


3507 


4257 


4187 






5/2 


-817 


4071 


3493 


4271 


4188 


3d^4p 


fijjo 


1/2 


-814 


4016 


3458 


4195 


4137 






3/2 


-815 


4037 


3467 


4267 


4156 






5/2 


-812 


3982 


3511 


4278 


4183 






7/2 


-816 


3989 


3451 


4280 


4129 






9/2 


-820 


4003 


3463 


4302 


4153 



III. RESULTS AND DISCUSSION 

CrII has five valence electrons, and these have a sig- 
nificant impact on the form of the basis orbitals. For the 
CI-I-MBPT method to work well it is important to have 
good initial orbitals, and so our Dirac-Fock and subse- 
quent B-spline codes include the Bd** or id^ configuration 
in the core, as described in Section Hi Al (Our calcula- 
tions show that saturation of the CI can be met satis- 
factorily in both V^~^ and V'^ approximations.) In the 
CI-f-MBPT code, the 3(i orbitals are then stripped from 
the core and become valence orbitals for the purposes of 
both the CI and MBPT components of the calculation. 
In this way excitations from the id shell are treated non- 
perturbatively. 

To calculate the MBPT diagrams one must include 
the change in effective core-potential, V^~^ — V'^ or 
V^~^ — V^~^ depending on the initial approximation. 
MBPT diagrams that include this interaction are known 
as subtraction diagrams, and in our calculation they are 
huge. There are three subtraction diagrams in S(i) and 
two in S*^^^ (these are shown in Figs. 2 and 4 of [24, re- 
spectively). When the subtraction diagrams in E^-*^^ are 
included, they significantly and adversely affect the en- 
ergies obtained, as can be seen in the CI-t-E^^-' columns 
(labelled "-t-E^i)") of Tables U and HH However, it turns 
out that these adverse effects are nearly completely com- 
pensated when E*^^-' and S*^'^^ are also included. When all 
second-order diagrams are included consistently the en- 
ergies and wavefunctions are improved by the addition. 

One might consider what happens if the subtraction 
diagrams are simply neglected from the calculation. In- 
deed the CH-S(i) energies are improved. However when 
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FIG. 1. Two large diagrams that affect the ground state 
multiplet that partially cancel. The labels refer to their des- 
ignation in (15i] . 



3d 3d 3d 3d 




2 (a) 3 (a) 

E^^) is added (either with or without the two-valence- 
electron subtraction diagrams) the energy levels obtained 
are again in very poor agreement with experiment. Thus 
it is not only inconsistent to leave out the subtraction 
diagrams, but it gives very poor results when all other 
second-order MBPT terms are included. 

The behaviour can be explained by examining the form 
of the S^i) and E^^) 

diagrams. For example, consider 
Fig. 2(a) and Fig. 3(a) from [23] with all external lines 
representing 3d electrons (see Fig. [1]). The subtraction 
diagram 2(a) has opposite sign to the zero multipole 
[k = 0) part of 3(a). This kind of cancellation is what 
finally suppresses the large subtraction diagrams, and is 
the reason why all second-order E diagrams must be in- 
cluded consistently. 

The energies obtained are slightly better in the V^^^ 
calculation. However, consider our calculations of the 
isotope shift constant k^ms , presented in Table IIIII We 
find that the pure CI calculations give different results 
in the V'^ and V'^~^ approximations, yet when the 
core-valence interactions are included self-consistently 
the agreement is much improved: the disagreement is re- 
duced from - 13% to less than 4% and in the Sd'^Ap 
transitions of astrophysical interest, more like '^1%. 

The calculation of g-values follow the same trend (Ta- 
ble |IV|. In this case the results are far less sensitive to 
details of the wavefunction: instead the relativistic ef- 
fects are determined by the form of the wavefunctions 
near the nucleus. We see very strong agreement between 
our and V^~^ results, especially after all E diagrams 
are included consistently. In the transitions of astrophys- 
ical interest, the different starting approximations leads 
to disagreements of the order ~ 25% in the pure CI case 
but - 5% when MBPT is included. 

Despite the consistency with respect to starting ap- 
proximation, our transition energies still differ from ex- 
periment. The most likely explanation is that we haven't 
taken full account of the relaxation of the core elec- 
trons. These have a strong effect on the 3d electrons 
via the exchange potential, yet relaxation of these or- 
bitals is only taken into account to second-order using 
perturbation theory. Ideally one would include them as 
valence electrons in the CI so that their relaxation could 



TABLE IV. g-factors and relativistic shifts, q (cm^^), for 
transitions to the ground state ^S^i2- Experimental g- 
factors are taken from Ref. [23]; calculated values are for the 
full CI-I-E method in the V'^~^ approximation. 

Level J ~g q (CI) q (CI-I-E) 









Expt. 


Calc. 


yN 


yN-l 


yN yN-l 


3d* 4s 




5/2 


1.669 


1.657 


-2483 


—2209 


O A OA 


—zooi 






7/2 


1.578 


1.587 


-2300 


—2034 


OOOQ 


— Zi40 




6po 


1/2 


—0.689 


—0.665 


—2052 


— 1748 


— iy / y 


— ioyo 






3/2 


1.124 


1.067 


-1959 


-1661 


— io ( 


— 1 / yz 






5/2 


1.314 


1.314 


-1807 


-1518 


— i / Uo 


— ioz4 






7/2 


1.378 


1.397 


-1597 


—1321 




1 qap: 

— ioyo 






9/2 


1.416 


1.434 


-1333 


-1073 




J- J-UU 






11/2 




1.454 


-1016 


-776 


-830 


-758 


3d*4p 


6po 


3/2 


2.382 


2.385 


-1607 


-1325 


-1489 


-1421 






5/2 


1.875 


1.880 


-1479 


-1209 


-1340 


-1280 






7/2 


1.710 


1.714 


-1281 


-1024 


-1117 


-1061 


3d*4p 


4po 


1/2 


2.844 


2.811 


-2146 


-1782 


-2122 


-2003 






3/2 


1.802 


1.786 


-1913 


-1517 


-1847 


-1704 






5/2 


1.624 


1.626 


-1089 


-1036 


-916 


-946 






1/2 


3.155 


3.186 


-1512 


-1316 


-1390 


-1357 






3/2 


1.824 


1.827 


-1373 


-1227 


-1232 


-1232 






5/2 


1.628 


1.634 


-1651 


-1204 


-1547 


-1396 






7/2 


1.577 


1.585 


-1417 


-1171 


-1281 


-1220 






9/2 


1.570 


1.552 


-1200 


-977 


-1036 


-990 



be treated non-perturbatively, however this is not pos- 
sible because the CI Hamiltonian size grows too large. 
In the V'^~^ approximation the core is more tightly 
bound, lessening the magnitude of relaxation terms. This 
likely explains the improved agreement with experiment. 

The 3p electrons also pose a potential problem for the 
isotope shift calculation. The scaled SMS operator that 
appears in ([3|) manifests itself in the dipole part (fc = 1) 
of the multipole expansion of Q [15J. Therefore it may 
be particularly affected if the 3p - 3d exchange terms are 
not adequately described by the method. It is for this 
reason that we conservatively use the difference between 
the CI and the CI+E calculation as an estimate of the 
uncertainty in /csmsj rather than the smaller difference 
between the and V^~^ calculations. 

We present our final values of relativistic shifts, for 
the transitions of astronomical interest (that is, those 
seen in quasar absorption spectra) in Table |Vl As ex- 
plained in Section mil we prefer our V^~^ approximation 
which gives better agreement with experimental transi- 
tion energies, however the difference between the and 
V'^~^ calculations is similar to the errors quoted. The 
uncertainty is estimated as the difference between our 
pure CI calculation and the full CI+MBPT calculation 
including all S diagrams. Our calculated g-values are 
seen to be in good agreement with the CI calculations 
of Rcf. ^30] . In Table IIVI we see that the experimental 
(jf-factors for these transitions are well reproduced by our 
calculation. There are no close levels in this case, so the 
methods of matching ^-factors ^3d\ are not required. The 
^^5/2 ^^'^ ^^5/2 transitions are strongly mixed (as can 
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TABLE V. Relativistic shifts, q, for transitions to ground 
state 3d^ 



Level J 


Ld \ q (cm 
(cm^""^) (A) this work 


Ref. [30] 


3d'*4p "P" 3/2 
5/2 
7/2 


48399 2066 -1421 (70) 
48491 2062 -1280(70) 
48632 2056 -1061 (70) 


1360(150) 
-1280(150) 
■1110(150) 


TABLE VL Calculated velocity structure in wavelength space 
of transitions to ground state in CrlL 


Upper level J 


A (A) (5A^'^VA) -c (kms-^) 

A = 50 A = 53 A = 5A 


3d''4p "P" 3/2 
5/2 
7/2 


2066 -0.535 (96) 0.252 (45) 
2062 -0.522 (84) 0.246 (39) 
2056 -0.535 (96) 0.252 (45) 


0.495 (89) 
0.484 (78) 
0.495 (89) 



be seen from the g-factors, which are well reproduced by 
our calculation), so if these transitions are ever seen in 
quasar absorption systems then a more careful analysis 
may be required. 

In a previous comparison between theory and exper- 
iment in SrII it was found that the SMS was underes- 
timated by theory at the ~ 30% level [3l|. This is a 
single- valence-electron ion and so there is no CI; rather, 
we added E^^' directly to our Dirac-Fock calculation. 
In the case of SrII the addition of S*^^-' did not ade- 
quately account for the effect of core-relaxation on the 
SMS. On the other hand, using the same method good 
agreement has been obtained between theory and exper- 
iment for Znll [l^, which we also treated as a single- 
valence-electron ion. Moreover, because the majority of 
the isotope shift in Cr II comes from the valence- valence 
contributions which are treated to all orders using CI, 
we have good reason to believe that our mass-shifts have 
been calculated with reasonable accuracy. 

We have also estimated the size of the field shift in 
these transitions using a small CI basis to estimate F and 
experimental values of (^(r^) taken from The field 

shift is expected to be small for a light element like CrII. 
Furthermore for the d — p transitions of astrophysical in- 
terest the orbitals do not overlap the nucleus strongly 
and so there is additional suppression. We find the field 
shift is of order ^ 0.010 kms^^ or smaller, which is much 
smaller than our uncertainty in /csms- We neglect it en- 
tirely. 

In Table IVII we present our isotope shift shift calcula- 
tions for astronomically relevant transitions of Cr II using 



the V^^^ results of Table UTTl in Again the uncer- 
tainty quoted is the difference between the pure CI and 
CI-I-MBPT calculations, i.e. the entire effect of S. This 
is very much larger than the difference between our V-^ 
and V^~^ calculations. We quote the velocity structure 
in wavelength space relative to the leading isotope ^^Cr. 
This is the preferred form for use in astronomy: the ve- 
locity shift is 5v = A"*'^^/A ■ c where A'^'^^ = A"* - A^^ 
and c is the speed of light in kms^^. 



IV. CONCLUSION 

We have shown that the CI-t-MBPT method can give 
transition energies in good agreement with experiment 
for low- lying transitions in CrII. Although the subtrac- 
tion diagrams are very large when the orbitals are cal- 
culated in the V^^^ or approximations, when all 
second-order MBPT diagrams are taken into account 
consistently the calculated energies are found to im- 
prove. This may help to direct future efforts using the 
CI-I-MBPT method in many-valence-electron ions such 
as Fell, which is of importance to studies of a- variation 
in quasar absorption systems. 

The SMS in CrII is found to dominate the isotope 
shift. For the %/2 M'^Ap transitions seen 
in quasar spectra they are five times the magnitude 
of the normal mass-shift and of opposite sign. The 
total mass-shift constant for these transitions, A^^ms ~t- 
^SMS — 3365 (112) GHz-amu (taking the J = 3/2 up- 
per level), is consistent with an earlier CI estimate of 
1900 (1200) GHz-amu [T?|, although clearly at the hmit of 
the uncertainty. The Cr II isotope shift is also quite large 
in comparison to many of the other isotope shifts used in 
the quasar analysis (although many are unknown). The 
velocity shift between the even isotopes 500 ms^^) is 
comparable to the isotope shifts of the AA2803 and 2796 
lines of Mgll (~ 850ms~"'^). Fortunately in the case of 
Cr II there are stable isotopes on either side of the lead- 
ing isotope, so one may hope that the total systematic 
shift due to variation of isotope abundances is small. 
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